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ABSTRACT 

We present a strong-lensing analysis of the galaxy cluster MS 1358.4+6245 (z = 
0.33), in deep 6-band ACS/HST imaging. In addition to the well-studied system at 
z = 4.92, our modelling method uncovers 19 new multiply-lensed images so that a 
total of 23 images and their redshifts are used to accurately constrain the inner mass 
distribution. We derive a relatively shallow inner mass profile, dlog S/dlog r ~ — 0.33± 
0.05 (r < 200 kpc), with a much higher magnification than estimated previously by 
models constrained only by the z — 4.92 system. Using these many new images we can 
apply a non-parametric adaptive-grid method, which also yields a shallow mass profile 
without prior assumptions, strengthening our conclusions. The total magnification of 
the z s — 4.92 galaxy is high, about a ~ lOOx over its four images, so that the inferred 
source size, luminosity and star- formation rate are about ~ 5x smaller than previous 
estimates, corresponding to a dwarf-sized galaxy of radius ~ 1 kpc. A detailed image 
of the interior morphology of the source is generated with a high effective resolution of 
only ~50 pc, thanks to the high magnification and to the declining angular diameter 
distance above z ~ 1.5 for the standard cosmology, so that this image apparently 
represents the best resolved object known at high redshift. 

Key words: dark matter, galaxies: clusters: individuals: MS 1358.4+6245, galaxies: 
clusters: general, galaxies: high- redshift, galaxies: formation, gravitational lensing 



1 INTRODUCTION 

The galaxy cluster MS 1358.4+6245 (also known as CL 
1358+62; hereafter MS 1358) was first identified by Zwicky 
& Herzog (1968), and later classified as an X-ray luminous, 
rich cluster by Luppino et al. (1991) based on observations 
from the Einstein Observatory Extended Medium- Sensitivity 
Survey (EMSS; Gioia et al. 1990, Stocke et al. 1991). This 
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cluster has been subject to intensive study, mainly due to its 
high X-ray luminosity and richness at an intermediate-high 
redshift (e.g., Henry et al. 1992, Carlberg et al. 1996, Pog- 
gianti & Barbaro 1996, Mushotzky & Scharf 1997, Nichol et 
al. 1997, Kelson et al. 1997), but became particularly famous 
when a pair of z — 4.92 lensed galaxies were uncovered in 
its inner region by Franx et al. (1997). Subsequently, other 
such high-z galaxies have been found in systematic surveys 
of galaxy clusters and in deep field imaging (e.g., Frye & 
Broadhurst 1998, Frye, Broadhurst & Benftez 2002, Kneib 
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et al. 2004, Stark et al. 2007, Bouwens et al. 2004, 2009b, 
2010, Bradley et al. 2008, Zheng et al. 2009). 

Magnified objects are particularly useful, being gener- 
ally bright enough to yield useful spectra, in particular the 
z — 4.92 system in MS 1358 which has an especially well- 
resolved large fold image straddling the tangential critical 
curve. Detailed spectroscopy by Franx et al. (1997; see also 
Swinbank et al. 2009) has revealed asymmetric and redward 
shifted Ly-alpha emission relative to metal absorption lines, 
which is claimed to imply this galaxy is suffering an outflow 
of metal enriched gas (Franx et al. 1997), similar to local 
dwarf starburst galaxies (Dekel & Silk 1986, Heckman, Ar- 
mus & Miley 1990). This spectral signature of gas outflow 
was shown to be a general property of z > 4 galaxies, in 
the larger sample of lensed galaxies behind eight massive 
clusters, by Frye, Broadhurst & Bemtez (2000), with impli- 
cations for the enrichment of the IGM in general (Scanna- 
pieco & Broadhurst 2001). 

In the past two decades other extensive studies have 
included or focused on MS 1358, ranging from X-ray to ra- 
dio measurements and scaling relations (e.g., Cagnoni, della 
Ceca & Maccacaro 1998, Cooray et al. 1998, Ettori & Fabian 
1999, Allen 2000, Wu 2000, Bohringer et al. 2000, Arabad- 
jis, Bautz & Garmire 2002, Brown et al. 2003, McCarthy et 
al. 2003, Egami et al. 2006, Laroque et al. 2006, Morandi, 
Ettori & Moscardini 2007), to intracluster content and evo- 
lution studies (e.g., van Dokkum et al. 1998, Kelson et al. 
2000a,b, 2006, Borgani et al. 1999, Fabricant, Franx & van 
Dokkum 2000, Ferreras & Silk 2000, Henry 2000, Kochanek 
et al. 2000, Tran et al. 2003, Postman et al. 2005, Holden 
et al. 2007), and arc statistics and searches (La Fevre et 
al. 1994, Hattory, Watanabe & Yamashita 1997, Luppino 
et al. 1999, Oguri, Lee & Suto 2003, Wambsganss, Bode & 
Ostriker 2004, Sand et al. 2005). 

Despite this broad study, only one multiply-lensed sys- 
tem has been hitherto identified in this cluster (the z = 4.92 
drop-out galaxy; Franx et al. 1997), so that no full SL anal- 
ysis has been possible given the wide degeneracy of models 
with so few images. Franx et al. (1997) have presented a SL 
mass model based on the z — 4.92 drop-out galaxy they un- 
covered, and various simplified mass models have been used 
by others (e.g., Allen 1998, Molikawa et al. 1999, Williams, 
Navarro & Bartelmann 1999), often based on isothermal po- 
tentials or spherical symmetry. Richard et al. (2008; see also 
Swinbank et al. 2009) have presented two z ~ 7.5 unverified 
multiply-lensed candidates for which Ly-alpha emission has 
not been detected. Moreover, we note that for a cluster at 
z c i — 0.33, the lensing-distance ratio for sources at z 3 = 4.92 
and z s ~ 7.5 differs by only ~ 2%, so that in order to mean- 
ingfully constrain the inner mass profile it is crucial to use 
other multiply-lensed systems at lower redshifts to expand 
the range of lensing-distances. 

Here we use our well-tested approach to modelling in or- 
der find a significant number of multiple images across the 
central field of MS 1358 so that the mass distribution and 
its profile can be well constrained. This method was devel- 
oped by Broadhurst et al. (2005a), and simplified further by 
Zitrin et al. (2009b) and has securely identified tens of mul- 
tiple images in high quality HST/ACS (Advanced Camera 
for Surveys) images, of background sources behind several 
clusters with deep ACS/HST imaging, including Abell 1689, 
C10024+17 and a sample of 12 MACS clusters at 2 > 0.5 



(Broadhurst et al. 2005a, Zitrin et al. 2009b, Zitrin et al. 
2010a). This is done with only six free parameters so that 
in practice the number of multiple images uncovered readily 
exceeds the number of free parameters as minimally required 
in order to obtain a reliable fit. 

This approach to lens-modelling is based on the rea- 
sonable assumption that mass approximately traces light. 
Recently we have independently tested this assumption in 
Abell 1703 (Zitrin et al. 2010b), by applying the non- 
parametric technique of Liesenborgs et al. (2006, 2007, 2009) 
for comparison, which is also employed here. This latter 
technique employs an adaptive grid inversion method and 
does not require prior assumptions regarding the mass dis- 
tribution, relying only on the images we have identified and 
their redshift estimates. For meaningful constraints, this ap- 
proach to modelling requires many sets of multiple images, 
and over a wide range of background redshifts per cluster. 
This condition is well met for A1703 (Zitrin et al. 2010b) 
and this model-independent method yields a very similar 
mass distribution to our parametric technique in the case 
of A1703, supporting the assumption that mass generally 
traces light. Independently, it has been found that SL meth- 
ods based on parametric modelling are accurate at the level 
of a few percent in determining the projected inner mass 
(Meneghetti et al. 2010). 

The results of this work will be further combined with 
a wide-range WL data from deep Subaru imaging to provide 
the cluster mass profile out to the virial radius and beyond 
(Umetsu et al., in preparation). To date only a few clusters 
have been reliably analysed by combining both weak and 
strong lensing for a full determination of the mass profile and 
a definitive comparison with theoretical predictions, (e.g., 
Gavazzi et al. 2003, Broadhurst et al. 2005b, 2008, Merten 
et al. 2009, Newman et al. 2009, Okabe et al. 2009 and ref- 
erences therein, Umetsu et al. 2010, Zitrin et al. 2010b). 
The upcoming multi-cycle HST program of cluster imaging 
(the CLASH program ) will provide a much more defini- 
tive derivation of mass profiles for a statistical sample of 
relaxed, X-ray selected clusters, combining high resolution 
space imaging with deep, wide-field ground based data, for 
a definitive determination of the equilibrium mass profiles 
of virilised clusters. 

Luppino et al. (1991) have measured the Brightest Clus- 
ter Galaxy (BCG) of MS 1358 at a redshift of z = 0.323. 
Other spectroscopic studies of many cluster member galax- 
ies suggest a similar cluster redshift of z ~ 0.33, which is the 
redshift adopted here (e.g., Fabricant, McClintock & Bautz 
1991, Carlberg et al. 1996, Yee et al. 1998, Fisher et al. 
1998). 

The paper is organised as follows: In §2 we describe the 
observations. In §3 we detail the SL analysis and in §4 we 
report and discuss the results, which are then summarised in 
§5. Throughout this paper we adopt a concordance ACDM 
cosmology with (Q m0 = 0.3, ^ao = 0.7, h = 0.7). With 
these parameters one arcsecond corresponds to a physical 
scale of 4.75 kpc for this cluster (at z = 0.33). The reference 
centre of our analysis is fixed at the centre of the BCG: RA 
= 13:59:50.55 Dec = +62:31:05.00 (J2000.0). 
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2 OBSERVATIONS AND PHOTOMETRIC 
REDSHIFTS 

MS 1358 was observed with the Wide Field Channel (WFC) 
of the ACS installed on HST in the framework of the ACS 
Guaranteed Time Observations (GTO; PI: Ford, H.; pro- 
gram IDs 9292, 9717, 10325) which includes deep observa- 
tions of several massive, intermediate-redshift galaxy clus- 
ters (Ford et al. 2003). Integration times of 7928, 5470, 
5482, 9196, 13552 and 17757 seconds were obtained through 
the F435W, F475W, F555W, F625W, F775W, and F850LP 
(Bg'Vr'i' z') filters, respectively. Some important aims of the 
GTO program were the determination of the mass distribu- 
tion of clusters for testing the standard cosmological model 
and to study distant, background lensed galaxies for which 
some of the very highest redshift galaxies are found because 
of high magnification by massive clusters. 

Mass models and detailed lensing analyses have been 
presented for most of the GTO clusters (e.g., Broadhurst 
et al. 2005a, Elfasdottir et al. 2007, Jee et al. 2007, Lemze 
et al. 2008, 2010, Limousin et al. 2008, Richard et al. 2009, 
Umetsu et al. 2010, Zitrin et al. 2009b, 2010b, Coe et al. 
2010, Medezinski et al. 2010). In this work, we present a 
strong lensing (SL) analysis of the ACS images of MS 1358, 
and detect several multiple-image systems so the mass pro- 
file of this cluster can be reliably constrained as well. 

The ACS images were initially reduced, processed, and 
analysed by APSIS, the ACS GTO pipeline (Blakeslee et 
al. 2003). An optimal \ 2 detection image was created as a 
weighted sum of all filters, each divided by its background 
RMS. Objects were detected and photometry obtained using 
SExtractor (Bertin & Arnouts 1996). 

The light of cluster galaxies was carefully modelled in 
each filter and subtracted from the images. This improves 
both the detection and photometry of lensed background 
objects. Based on this Bg'Vr'i' ' z photometry, we obtain 
photometric redshifts using BPZ (Bemtez 2000, Bemtez et 
al. 2004, Coe et al. 2006). The distances to the galaxies are, 
of course, key ingredients to the lens model. 



3 STRONG LENSING MODELLING AND 
ANALYSIS 

We apply our well tested approach to lens modelling, which 
has previously uncovered large numbers of multiply-lensed 
galaxies in ACS images of Abell 1689, C10024, and 12 high-z 
MACS clusters (respectively, Broadhurst et al. 2005a, Zitrin 
et al. 2009b, Zitrin & Broadhurst 2009, Zitrin et al. 2009a, 
2010a, b). Briefly, the basic assumption adopted is that mass 
approximately traces light, so that the photometry of the red 
cluster member galaxies is used as the starting point for our 
model. Cluster member galaxies are identified as lying close 
to the cluster sequence by the photometry described in Sj2] 
We approximate the large scale distribution of cluster 
mass by assigning a power-law mass profile to each galaxy 
(see Figure [T]), the sum of which is then smoothed (see Fig- 
ure The degree of smoothing (S) and the index of the 
power-law (q) are the most important free parameters de- 
termining the mass profile. A worthwhile improvement in 
fitting the location of the lensed images is generally found 
by expanding to first order the gravitational potential of 
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Figure 1. The starting point of the mass model, where we de- 
fine the surface mass distribution based on the cluster member 
galaxies (see ^3|. Axes are in ACS pixels (0.05" /pixel). 
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Figure 2. The resulting smooth mass component of the mass 
model (see {J3jl. Axes are in ACS pixels (0.05" /pixel). 

this smooth component, equivalent to a coherent shear de- 
scribing the overall matter ellipticity, where the direction of 
the shear and its amplitude are free parameters, allowing 
for some flexibility in the relation between the distribution 
of DM and the distribution of galaxies, which cannot be ex- 
pected to trace each other in detail. The total deflection field 
St(0), consists of the galaxy component, a ga i(8), scaled by 
a factor K ga i, the cluster DM component 5dm (0), scaled by 
(1-Kg a i), and the external shear component a ex (0): 

a T (6) = Kgaitigaltf) + (1 - K gal )a DM (0) + a ex (6), (1) 

where the deflection field at position 6 m due to the external 
shear, d ex {0 m ) = (a ex ,x,aex, v ), is given by: 

ctex,x{0 m ) = \l\ cos(2(j) 1 )Ax m + |7| sin(20 7 )Ay m , (2) 

OLex, y {0 m ) = I7I sin(2<?!) 7 )Aa; m — | — ^ | cos(2^ 7 ) Ay m , (3) 

where (Ax m , Ay m ) is the displacement vector of the po- 
sition 6 m with respect to a fiducial reference position, which 
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we take as the lower-left pixel position (1, 1), and <j> y is the 
position angle of the spin-2 external gravitational shear mea- 
sured anti-clockwise from the a;-axis. The normalisation of 
the model and the relative scaling of the smooth DM compo- 
nent versus the galaxy contribution brings the total number 
of free parameters in the model to 6. This approach to SL 
is sufficient to accurately predict the locations and internal 
structure of multiple images, since in practice the number 
of multiple images uncovered readily exceeds the number of 
free parameters thus fully constraining them. 

In addition, two of the six free parameters can be pri- 
marily set to reasonable values so only 4 of these parame- 
ters have to be constrained initially, which sets a very re- 
liable starting-point using obvious systems. The mass dis- 
tribution is therefore primarily well constrained, uncovering 
many multiple-images which can then be iteratively incor- 
porated into the model, by using their redshift estimation 
and location in the image-plane. 

Firstly we use this preliminary model to lens the more 
obvious lensed galaxies back to the source plane by subtract- 
ing the derived deflection field, and then relens the source 
plane to predict the detailed appearance and location of ad- 
ditional counter images, which may then be identified in 
the data by morphology, internal structure and colour. We 
stress that multiple images found this way must be accu- 
rately reproduced by our model and are not simply eyeball 
"candidates" requiring redshift verification. The best fit is 
assessed by the minimum RMS uncertainty in the image 
plane: 

RMSf mages = / J((a;< - £i) 2 + {Vi - Vi) 2 ) I N image s, (4) 

i 

where x i and y i are the locations given by the model, and 
Xi and yi are the real image locations, and the sum is over 
all Ni ma g ea images. The best-fit solution is unique in this 
context, and the model uncertainty is determined by the 
location of predicted images in the image plane. Impor- 
tantly, this image-plane minimisation does not suffer from 
the well known bias involved with source plane minimisa- 
tion, where solutions are biased by minimal scatter towards 
shallow mass profiles with correspondingly higher magnifi- 
cation. 

The model is successively refined as additional sets of 
multiple images are incorporated to improve the fit, impor- 
tantly using also their redshift information for better con- 
straining the mass slope through the cosmological relation 
of the di s /d s growth. 

In order to examine our results, we also apply the 
adaptive-grid, non-parametric modelling method of Liesen- 
borgs et al. (2006, 2007, 2009), which makes no prior as- 
sumptions regarding the mass distribution. In a recent work, 
we have made a more thorough comparison of these meth- 
ods by analysing the well studied cluster Abell 1703 (Zitrin 
et al. 2010b). We showed that the results of these two dis- 
tinct approaches are very similar when sufficient data con- 
straints are available. The effective area modelled with the 
non-parametric technique is limited to roughly the critical 
area, beyond which there are very few multiple-images that 
can be used to reliably constrain the fit. Some discrepancy 
may be expected with regards to the mass distribution be- 
tween the two methods, especially when not enough con- 
straints are available. However, aside from minor local trade- 



offs originating in the mass-sheet degeneracy or in the inclu- 
sion of cluster members, the overall profile can be expected 
to be similar in order to maintain the cosmological lensing- 
distance ratio of the various systems. 



4 RESULTS AND DISCUSSION 

4.1 Multiple-Images, Mass Model and Critical 
Curves 

In addition to the four previously-known lensed images of 
the z 3 = 4.92 dropout galaxy (Franx et al. 1997), our mod- 
elling technique has uncovered 19 new multiply-lensed im- 
ages in the central field of MS 1358, belonging to 7 new 
systems. We have made use of the location and photometric 
redshift information of these images to fully constrain the 
mass model. We find that the critical curves for a source 
at z a = 4.92 (system 1) enclose an area with an effective 
Einstein radius of te — 21 ± 3", or ~100 kpc at the red- 
shift of the cluster. This critical curve we find (see Figure 
[3]) encloses a projected mass of M ~ 6.1 ± 0.8 x 1O 13 M . 
For a source redshift of z s ~ 2 the effective Einstein ra- 
dius is te = 13 ± 2", enclosing a projected mass of M ~ 
2.7 ± 0.2 x 10 13 M Q around the BCG. For general compari- 
son, this is in very good agreement with the Einstein radius- 
mass scaling relation for a source at z s — 2, found in Zitrin 
et al. 2010a (taking into account also the different lens dis- 
tances; see Figure 27 therein). The corresponding critical 
curves are plotted on the cluster image in Figure [3] along 
with the multiply-lensed systems. The resulting mass distri- 
bution and its profile are shown in Figures |4]and[5] 

Note that another small critical curve is formed around 
the group of bright cluster galaxies ~ 1' to the south, whose 
prominence is hard to determine due to lack of multiple- 
images at that region. We note that if this southern clump 
is somewhat more massive than presented here, the critical 
curves may merge with the central main critical curves to 
form a larger more elongated critical region. This however 
does not seem probable as only few arcs are seen in between, 
and thus rule out the existence of such an extended critical 
curve. 

Several mass models have been created in earlier work, 
describing the critical curves and mass distribution based 
only on system 1 (z = 4.92). Franx et al. (1997) used isother- 
mal potentials with which they were able to reconstruct the 
images of this system. They find an Einstein radius of 21", 
similar to our result. Allen (1998) has calculated a total 
mass of 8.27 x 1O 13 M0 with 20% uncertainty, within a ra- 
dius of 121 kpc (25.5"), in good agreement with our result 
considering the Einstein radius difference. Other mass mod- 
els made for this cluster based on projected WL profiles and 
circular symmetry, derive similar or slightly smaller values 
(see Diaferio, Geller & Rines 2005, e.g., Hoekstra et al. 1998, 
2007, Arabadjis, Bautz & Garmire 2002, Takahashi & Chiba 
2007). 

It should be stressed that the multiple-images found 
here are accurately reproduced by our model and are not 
simple identifications by eye. The parametric method of 
Zitrin et al. (2009b) has been shown in many cases to have 
the predictive power to find multiple images in clusters. Due 
to the small number of parameters this model is initially 
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Figure 3. Galaxy cluster MS 1358 (z = 0.33) imaged with Hubble/ACS. We number the multiply-lensed images used and uncovered 
in this work. The overlaid white critical curve corresponds to system 1, the red drop-out galaxy at z B = 4.92 (found by Franx et al. 
1997), enclosing a critical area of an effective Einstein radius of ~ 100 kpc at the redshift of this cluster. Interior to this, the blue critical 
curves correspond to systems 3,4,5, at a lower redshift of z a ~ 2. Additional less secure candidate systems (marked with "c") have similar 
appearance but are not well reproduced by our best fit solution, and were not used to constrain the model. We also show the approximate 
location of two unverified candidate z s ~ 7.5 multiple-systems, marked with "z" (Richard et al. 2008). These were not used to constrain 
the model but their locations and claimed high redshift are plausible in the context of our model, see §4.1| 
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Table 1. Multiple-image systems and candidates used and uncovered by our model. For more detailed information on each system 
and other optional members see the corresponding subsection. The columns are: arc ID; RA and DEC in J2000.0; 95% confidence level 
minimal and maximal photo-z; best photo-z; z mo< j e !j the redshift predicted by the mass model; comments. Note that system 1 was 
uncovered by Franx et al. (1997) who measured their redshift spectroscopically, which is the value given below. Note also that unusually 
large errors in the photo-z imply a bimodal distribution. In such cases the values which agree with the SL model can be different than 
specified in the best photo-z column, as they arise from another peak in the distribution. Such cases are specified in the comments. 
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well-constrained enabling a reliable identification of other 
multiple-images in the field, which are then used to fine- 
tune the mass model. We now detail each multiply-lensed 
system in turn, as listed in Table [l] 



System 1 : A high-redshift dropout galaxy at z s = 4.92. 
The four images of this multiply-lensed source were found 
by Franx et al. (1997) who also measured its redshift. Two of 
the multiple-images of this system (1.1/1.2), form a promi- 
nent red fold-arc about ~ 21" south-west of the BCG, next 
to another bright cluster member (see Figure [3|. Image 1.3 
is a smaller image on the other side of that cluster member, 
and image 1.4 is ~ 29" east of the BCG. These images are 
well reproduced by our model, as seen in Figure [6] Accord- 
ing to our model, the angular area of the source is ~ O.lfj" 
(see Figure FfJ, and is therefore magnified about a 100 times 



in area, summed over all four images which in total subtend 
~ 10D" in the image-plane (see also E 4.2 I . 



System 2 : A long faint arc ~ 21" north of the BCG, 
next to a bright cluster member. Our model reproduces this 
arc accurately at a redshift of z s ~ 3.4, similar to the pho- 
tometric redshift of image 2.1. The photometric redshift of 
image 2.2 is double peaked, with peaks around 0.5 and 3.4, 
the latter in agreement with our model. In addition, the 
clearly lensed arc lies away from the BCG and should have a 
relative high redshift in order to be lensed, further strength- 
ening our match. Another similar looking candidate is seen 
to the right of image "clc" (see Figure [3JI, but as can be 
seen the critical curves do not pass between these images 
and thus this third candidate is not probable in the context 
of our model. 

Systems 3, 4, and 5 : These three systems consist of typ- 
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ical blue and white arclets seen in large amounts in lensing 
clusters. Systems 3-5 follow a similar symmetry, as can been 
seen in Figure [3] One image of each system appears ~ 10" 
north of the BCG, the second image appears ~ 15" east of 
the BCG and the third image ~ 20" to its west. These three 
systems have similar typical photometric redshifts of 2 ~ 2 
(see Table[l]|, and due to their close vicinity and similar ap- 
pearance were hard to match accurately, especially in the 
case of the third images to the west. We therefore note that 
there can be an irrelevant mix-up with regards to the third 
image of each of these three systems, but since the distances 
between them are small this does not affect the resulting 
mass model. In addition, due to a minor, northwards off- 
set in the image-plane reproduction of images 4.3 and 5.3, 
we acknowledge the possibility that there are other similar- 
looking spots which may correspond to these systems, but 
have a negligible effect on the mass model. 

System 6 : Two blue faint images with mirror symmetry 
lying on two sides of the critical curve for a source at 2 ~ 1.8 
by our model. Their photometric redshifts are ~ 2.0 ± 0.4 
and ~ 1.77lg'4, respectively, where the first may be affected 
by a nearby galaxy light. Our model predicts the extra small 
and faint spec on the other side of the cluster about ~ 20" 
south-west of the BCG, as identified in the data. 

System 7 : Two greenish small images on two sides 
of a bright galaxy. These lack photometric redshifts due to 
the vicinity of the galaxy, and our model reproduces them 
accurately at 2 ~ 1.8. 

System 8 : Three blue and white relatively bright im- 
ages of a 2 ~ 2.1 galaxy, following similar symmetry as sys- 
tem 1. Our model reproduces these images very well at this 
redshift, as seen in Figure [8] Note in Figure [I] the highly 
magnified region near the location of image 8.1, which ex- 
plains its larger appearance relatively to images 8.2 and 8.3. 

Other candidate systems and image identification un- 
certainty: A similar symmetry to that systems 2-5 follow, 
can be expected in the northern region of the cluster be- 
tween the BCG and system 2. Moreover, several images in 
the west side of the cluster (marked in Figure [T]) resem- 
ble their corresponding image candidate near system 2, and 
have reasonably similar photo-2's. However, according to our 
model these candidate systems are not being lensed in full 
or in their expected photometric redshifts and we therefore 
exclude them from our secure identification. 

System "cl" consists of two bright blue-white arclets 
which seem to have mirror symmetry ("cla" and "clb"), 
and an additional similar looking arclet, "clc", ~ 10" south- 
west of them. These have a photometric estimate of z 3 ~ 2 
and as can be seen in Figure [l] are north of the critical 
curve for a such a source, so that they are not lensed by 
our model. This may be artificially overcome if one boosts 
the mass obtained by the nearby galaxy, but then system 2 
would not correspond to its z s ~ 3.4 photometric redshift. In 
addition, both this system, and system 2, have candidates 
to the west, which are not generated by our model, since 
there is no extra apparent mass nor galaxies between these 
locations, so the critical curve does not pass between them. 
One can significantly boost the nearby galaxies to create 
a mass lump which might produce images at that region; 
however, it should be noted that such actions are rarely 
justified and have not been critically required in any of the 
many clusters we have analysed to date. We conclude that 
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Figure 4. 2D surface mass distribution (re), in units of the critical 
density (for z B = 4.92), of MS 1358. Contours are shown in linear 
units, derived from our mass model constrained using the many 
sets of multiply-lensed images seen in Figure [3] Axes are in ACS 
pixels (0.05" /pixel), and a 20" bar is overplotted. 

despite the similar appearance this system is not probable 
in the context of our model. 

System "c2" is a much more probable candidate con- 
sisting of four red drop-out arclets. Our model reproduces 
and therefore identifies arcs "c2c" and "c2d" as the same 
system in any redshift in the range 3 to 5. Arcs "c2a" and 
"c2b" are obtained as part of this system only in redshifts 
higher than ~ 4.5, corresponding to the photometric red- 
shift of arcs "c2a" and "c2b", of ~ 4.5. Still, we decided not 
to include it as a secure system since the redshift agreement 
is marginal and the distance in the image-plane reproduc- 
tion of arcs "c2a" and "c2b" (from their observed location) 
is larger than our typical uncertainty. 

Systems "zl" and "z2" are multiple-images of z 3 ~ 7.5 
dropout candidates, claimed as possibilities by Richard et 
al. (2008). In Figure [3] we mark their approximate locations 
generated by our model (see Figure 12 in Richard et al. 
2008). We did not use these candidate images to constrain 
the model, but verified that they are plausible high-2 objects 
in the context of our model (for z c i = 0.33 the lensing- 
distance of sources at z a ~ 7.5 is only ~ 2% bigger than the 
lensing distance of the z s = 4.92 galaxy, so that redshifts 
throughout that range are plausible in that context). 

4.2 Mass Profile and High Magnification 

We mentioned in the preceding sections that the profile can 
only be accurately constrained by incorporating the cosmo- 
logical redshift-distance relation, i.e., the lensing distance of 
each system based on the measured spectroscopic or photo- 
metric redshifts. In so doing we normalise our mass model 
to system 1, so that the normalised scaling factor, f(di a /d a ), 
is equal to 1. We then make use of the lower-2 systems, in 
particular the 2 ~ 2 systems (number 3,4,5 and 8) and the 
2 ~ 3.4 system (number 2), whose photometric redshifts we 
find most credible to constrain the profile. 

We examine how well the cosmological relation is repro- 
duced by our model, accounting for all systems with photo- 
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Figure 5. Radial surface mass density (re) profile in units of the 
critical surface density (for a fiducial redshift of z s = 4.92), de- 
rived using the sets of multiple images shown in Figure [3] The 
cyan shaded curve is derived by our modelling method (Zitrin 
et al. 2009b), while the dashed black lines are derived using 
the non-parametric technique of Liesenborgs et al. (2006). As 
can be seen, the profiles agree very well in the region of in- 
terest, and are rather flat. We measure a profile slope of only 
d log S/d log 9 ~ —0.19 ± 0.05 within the effective Einstein radius 
for system 1 (r^ = 21 ± 3; ~ 100 kpc; marked with a vertical 
black line), increasing to d log T,/d log 9 ~ —0.33 ± 0.05 within 
twice the critical range, ~ 200 kpc. Note that the inner profiles 
are at a level close to the critical density (±0.25) for system 1. 
For comparison, overplotted with blue dash-dotted lines are kappa 
profiles for three other well-known clusters: Abell 1703, CI 0024, 
and Abell 1689, after correcting for the different lens distances. 
These were slightly shifted upwards on the y-axis to allow a clean 
view and to better demonstrate the slope difference of MS 1358 
from these typical lensing clusters which exhibit a common mass 
slope of dlog £/dlog# ~ —0.5. 
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Figure 6. Reproduction of system 1 by our model, by delensing 
image 1.1 into the source plane, and then relensing the source 
plane pixels onto the image plane. Our model clearly reproduces 
very accurately the second half of the arc, 1.2, and the internal 
structure of other images in this system including the prominent 
internal HII regions within the source. 




Figure 7. Reproduction of the z s = 4.92 source galaxy by our 
model, using the RIZ images. The image is reproduced by de- 
lensing arc 1.1/1.2 into a high-resolution source-plane chosen to 
accommodate the high magnification at its location. Two blue 
lower-redshift objects got delensed to the source plane in the pro- 
cedure. Note that the source shape is very similar to that found 
in Swinbank et al. (2009), though due to the magnification differ- 
ence it is smaller by a factor of about 4, and seem less internally 
stretched. The color-coding was slightly modified to obtain a bet- 
ter view of the internal details, and some noise was removed from 
the edges of the image. 
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Figure 8. Reproduction of system 8 by our model, by delensing 
image 8.1 into the source plane, and then relensing the source 
plane pixels onto the image plane. Our model clearly reproduces 
very accurately the other images in this system. The clear size dif- 
ference of image 8.1 compared with images 8.2 and 8.3 (zoomed-in 
here for better view) is another indication of the high local mag- 
nification in that region. 

z's, as shown in Figure [9] Clearly the redshifts of these sys- 
tems verify very well that the predicted deflection of the best 
fitting model at the redshift of each of these systems, lies 
accurately along the expected cosmological relation, with 
a mean deviation of only A/ < 0.01 (see Figure [9j|, and 
X 2 — 0.1 for the best model, considerably strengthening the 
plausibility of our approach to modelling in general. 

We find that unlike other well-known lensing clusters 
such as Abell 1689, Abell 1703 or C10024+17, which all have 



Strong- Leasing Analysis of MS 1358.4 +624 5 9 



a typical profile slope of dlogE/dlog# ~ —0.5 (see Broad- 
hurst et al. 2005a, Zitrin et al. 2009b, 2010b), MS 1358 has 
a relatively shallow profile slope of only d log E jd log 8 ~ 
— 0.19±0.05 in its main critical region of ~[3",2l"], increas- 
ing to dlog£/a!log# ~ —0.33 ± 0.05 in twice the critical 
region ~[3",42"] (~ 200 kpc), as can be seen in Figure 
[5] We compare our profile with the profile derived by the 
assumption-free method of Liesenborgs et al. (2006), over- 
plotted as well in Figure [5] As can be seen the profiles are 
similarly shallow and in good agreement with each other. 
This may be expected, since the profile must be adjusted to 
maintain the lensing-distance ratio for each lensed system, 
manifesting the importance of the images found in this work 
for constraining the mass profile of this cluster. 

The shallow profile results in high magnification, as 
kappa is close to unity over a wide central area and thus the 
magnification is boosted. This phenomenon we have revealed 
recently in another even shallower profile cluster, MACS 
1149.5+2223 (Zitrin & Broadhurst 2009) which hugely mag- 
nifies a distant spiral galaxy by a factor of ~ 200 x creat- 
ing the largest lensed images known, as the slope is even 
shallower in this case related to the unrelaxed nature of 
that cluster. Another work based on spectroscopic red- 
shifts claims a somewhat more moderate slope of up to 
dlogE/dlogfl ~ -0.3 for MACS 1149.5+2223 (Smith et 
al. 2009) but still significantly lower than typically found 
in other rich lensing clusters. In MS 1358 the high mag- 
nification is inferred by the distribution of multiple-images 
and their redshifts, and can account for the large bright im- 
ages of system 1, at z s = 4.92 (see Figure [3|. We find that 
the source galaxy is magnified ~ 80 x into arcs 1.1/1.2, and 
about 5 x — 15x into each of arcs 1.3 and 1.4, so that the 
source is magnified in total about a 100 times, and is there- 
fore one of the most highly-magnified distant objects known. 

The magnification values are hard to determine pre- 
cisely since the magnification is a very sensitive function of 
the local mass gradient so that in practice shallower models 
yield much higher magnifications. For example, Franx et al. 
(1997), have used isothermal potentials to find a magnifica- 
tion of x 5 — x 11 at the eastern part of the bigger arc (image 

1.2 here), and Swinbank et al. (2009; see also Richard et al. 
2008) have found a revised "luminosity-weighted" magnifi- 
cation of ~ 12.5, while we find in that region an average 
magnification of ~ x25 — x50. It is well understood now by 
most practitioners that the mass profile cannot be appro- 
priately constrained when too few multiple-images are used, 
or when they span only a narrow redshift range. Here, our 
identification of many multiple-images spread over a wide 
range of redshifts allows for the first time an accurate de- 
termination of the profile slope of MS 1358, revealing the 
rather shallow central mass profile. 

4.3 Luminosity of the z = 4.92 Source 

Due to the large difference in magnification with models 
published earlier (Franx et al. 1997, Swinbank et al. 2009), 
and due to the importance of the object at z s = 4.92, we cal- 
culate our corrections to the original z = 4.92 source lumi- 
nosity. Firstly, we simply divide the bolometric luminosity 
published by Franx et al., ~ 3 x 10 11 Lq, by the approxi- 
mate magnification ratio of the two mass models, to obtain 
a corrected source bolometric luminosity of ~ 6 x 10 10 I/q. 
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Figure 9. Growth of the scaling factor f(di a /d a ) as a function 
of redshift, normalised so f=l at z = 4.92. Plotted lines are the 
expected ratio from the chosen specified cosmological model. The 
circles correspond to the multiple-image systems reproduced by 
the parametric mass model, versus their real photometric redshift. 
The data follow very well the relation predicted by the standard 
cosmological model (mean deviation of only Af < 0.01, and x 2 — 
0.1 for this fit). We take an average photo-z for each system, where 
there is an estimation difference between the same images of the 
same system. 

Independently, we find that the source AB magnitude is 
1814, ab = 25.67 and I775,ab = 25.76 (after accounting for 
the magnification), significantly fainter than lau.AB = 24.0 
or I775,as = 24.94, as published by Franx et al. and Swin- 
bank et al., respectively. We then use eq. 2 of Bouwens et 
al. (2009a) to convert the corrected bolometric luminosity 
into an absolute magnitude in the rest-frame UV. This re- 
lation is known to hold from z ~ 2.5 to at least z ~ 4, 
and is expected to be approximately valid also at z = 4.92. 
Following our model, the absolute magnitude of the source 
is My VAB ~ —19.5, ~ 4 times a more common value for 
objects at z ~ 5 according to up-to-date high-z luminosity 
functions (e.g., Bouwens et al. 2007, van der Burg, Hilde- 
brandt & Erben 2010). 

4.4 Mass- To-Light Ratio 

We calculate here the central Mass-To-Light ratio using 
the Luminous Red Galaxies (LRG) template described in 
Bem'tez et al. (2009) to convert fluxes into luminosities, and 
measure an M/ Lb ratio of ~ 160 ± 30 (M/L)q within the 
critical curves for z a = 4.92 (see Figure 

!}. We compare 

this value to the M/L versus Einstein radius relation found 
in Zitrin et al. (2010a) for a sample of 12 clusters. After 
correcting for the difference in the lens redshift, this value 
tightly follows the relation they presented (see Figure 28 
therein), which predicts a similar ratio. 

We may also examine in detail the M/L of the lu- 
minous elliptical galaxy which splits the images of the 
z s — 4.92 galaxy arcs 1.1/1.2 and 1.3 (visible in Figure 
[3] RA=13:59:49.47 DEC= +62:30:47.55). According to our 
model this cluster member has a mass of 1.451q'4q x 1O 11 M0 
within a radius of 3" (~ 14 kpc), after subtracting the clus- 
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ter smooth DM component. We compare this result with the 
mass inferred by the virial theorem through the velocity dis- 
persion, and in the same aperture. For example, Kelson et al. 
(2000a, b, 2006) measured a velocity dispersion of a ~ 220 
km/s for this galaxy, which corresponds to ~ 1.6 x 1O 11 M0 
within 3", in fair agreement with our result (we note however 
that a smaller galaxy may be hidden in this bright elliptical 
light, slightly affecting the comparison). The light enclosed 
within 3" is ~ 6 x 10 10 Lq, which yields a small M/Lb ratio 
of ~ 3 (M/L)q, close to the purely stellar value indicat- 
ing no significant DM out to ~ 14 kpc. The effective radius 
of this galaxy was found to be ~ 1" (4.75 kpc; Kelson et 
al. 2000a,b, 2006) so that in practice we measure in a ra- 
dius bigger three times than the effective radius, where the 
DM should be dominant. However, a lower M/Lb ratio is 
not surprising, because we subtracted the cluster DM com- 
ponent from our mass measurement, and the halo of this 
galaxy should already be stripped by the significant tidal 
forces in the central cluster region (if not earlier during the 
hierarchical evolution of the cluster through mergers of sub- 
clumps). It is known that about 70% of the DM is stripped 
from central galaxies also in simulations of massive clusters 
(e.g., Nagai & Kravtsov 2005). 



5 SUMMARY 

In this work we have presented a detailed lensing analysis 
of the galaxy cluster MS 1358 in HST/ACS images. Our 
well-established modelling method (Broadhurst et al. 2005a, 
Zitrin et al. 2009b, Zitrin et al. 2010a, b) has revealed the 
rather shallow mass distribution of the central region by 
uncovering 19 multiply-lensed images which eluded previous 
detection, so that in total 23 images of 8 different sources 
were used to fully constrain the fit. Though more lensed 
candidates might be found in this field with further careful 
effort, our minimalistic approach to lensing involves only 6 
free parameters so that the resulting model is clearly fully 
constrained by these multiple images. In addition we have 
uncovered several other lensed candidates, for which it would 
be interesting to obtain spectroscopic redshifts (as well as to 
the systems identification presented here) in order to further 
establish the results of this work. 

The photometric redshifts of the newly found arcs 
throughout the central region enable the determination of 
the inner mass profile of MS 1358, through the cosmological 
lensing-distance ratio, and imply a shallow mass distribution 
manifested also in our modelling iterations so that mod- 
els spanning almost the full parameter space yield similar 
shallow mass profiles with typical slopes of d log E jd log 9 ~ 
-0.3±0.1 (r < 200 kpc). We further tested our results with 
a non-parametric adaptive-grid method (Liesenborgs et al. 
2006, 2007, 2009), which yields a similarly shallow profile 
hence strengthening the conclusions of this work. A shallow 
mass distribution translates into a boosted magnification in 
the central region so that the z — 4.92 source galaxy is mag- 
nified in total ~ 100 times, and is therefore one of the more 
highly-magnified, distant objects known. 

Our magnification values for this system are higher by a 
factor of ~ 5 than previous estimates (e.g., Franx et al. 1997, 
Swinbank et al. 2009) based on models constrained only by 
the z = 4.92 system. Hence the source size, bolometric lu- 



minosity, and star- formation rate are correspondingly ~ 5x 
smaller, more typical of faint field galaxies based on current 
high- 2 luminosity functions (e.g., Bouwens et al. 2007, van 
der Burg, Hildebrandt & Erben 2010), and with a physical 
scale of a dwarf galaxy of radius ~ 1 kpc. A detailed im- 
age of the interior morphology of the source is generated by 
delensing the most magnified image (1.1/1.2), resulting in 
a high effective resolution of only 50 pc, thanks to the high 
magnification and to the declining angular diameter distance 
above z ~ 1.5 for the standard cosmology, so that this im- 
age apparently represents the best resolved object known at 
high redshift. 
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